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ABSTRACT 

An ideal coronagraph with a band-limited image mask can efficiently image 
off-axis sources while removing identically all of the light from an on-axis source. 
However, strict mask construction tolerances limit the utility of this technique 
for directly imaging extrasolar terrestrial planets. We present a variation on the 
basic band-limited mask design — a family of "notch filter" masks — that miti- 
gates this problem. These robust and trivially achromatic masks can be easily 
manufactured by cutting holes in opaque material. 

Subject headings: astrobiology — circumstellar matter — instrumentation: adap- 
tive optics — planetary systems 



1. INTRODUCTION 

Direct optical imaging of nearby stars has emerged as a potentially viable method for 
detecting extrasolar terrestrial planets, buoyed by new techniques for controlling diffracted 
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and scattered light in high-dynamic-range space telescopes (see, e.g., the review by Kuchner 
& Spergel 2003). These techniques boost a telescope's ability to separate a planet's light 
from the light of its host star. At optical wavelengths, the Sun outshines the Earth by a 
factor of nearly 10^°; this contrast ratio is ~ 10^ times larger than the contrast ratio in 
the mid-infrared (Beichman et al. 1999; Des Marais et al. 2001). But to offset the higher 
dynamic range requirements of visible-light planet finding, optical techniques offer freedom 
from large, multiple-telescope arrays (Woolf 2003), cryogenic optics, and background light 
from zodiacal and exozodiacal dust (Kuchner & Brown 2000), while providing access to O2 
and O3 biomarkers (Traub & Jucks 2001; Des Marais et al. 2001), surface features (Ford, 
Seager & Turner 2001), the total atmospheric column density (Traub 2003), and even 
potentially the "red edge" signal from terrestrial vegetation (Woolf et al. 2002). 

Of the obstacles to achieving the necessary dynamic range in a single-dish optical tele- 
scope, the diffracted light background appears relatively manageable. For example, main- 
taining the scattered light background at the level of the expected signal from the planet 
poses a greater challenge; this task requires a r.m.s. wavefront accuracy of < 1 A (Kuchner & 
Traub 2002; Trauger et al. 2002a) over the critical spatial frequencies. However, techniques 
for managing the diffracted light may dictate the general design of a planet-finding telescope 
and the planet search and characterization strategy. 

Optical techniques for controlling diffracted light in planet-imaging telescopes have cen- 
tered on two main designs: specially shaped and/or apodizcd pupils (Spergel 2001; Nisenson 
& Papaliohos 2001; Kasdin, Spergel & Littman 2001; Dcbcs et al. 2002; Kasdin ct al. 2003) 
and classical coronagraphs (Lyot 1939; Nakajima 1994; Stahl & Sandler 1995; Malbet et al. 
1995; Kuchner & Traub 2002). Shaped and apodized pupils produce a point spread function 
whose diffraction wings are suppressed in some regions of the image plane. A classical coro- 
nagraph explicitly removes the on-axis light from the optical train by reflecting or absorbing 
most of it with an image mask and diffracting the remainder onto an opaque Lyot stop. 

Recently, Kuchner & Traub (2002) showed that a classical coronagraph performs best 
with a "band-limited" image mask. Different band-limited masks offer high performance for 
planet searching or planet characterization. For planet characterization, the sin^ amplitude 
transmissivity mask (sin'* intensity transmissivity) introduced in Kuchner & Traub (2002) 
can achieve 80% throughput for a planet at AX/D. With this high throughput, a 10 m by 
4 m telescope can detect a planetary biomarker in ~l/3 of the time needed by alternative 
designs (e.g., an 8 m square apodizcd aperture). A band-limited mask of the form 1— sine (see 
Table 1) has both excellent throughput and large search area. With any band- limited mask, 
an ideal coronagraph eliminates identically all of the on-axis light, though pointing errors 
and the stellar size contribute to a flnite leakage (Kuchner & Traub 2002). A band-hmited 
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mask can operate with a pupil of any shape as long as it has uniform transmissivity. 

But because they interact with focused starlight, all coronagraphic image masks face 
severe construction tolerances. Errors in the mask intensity transmissivity of ~ 10~^ on 
scales oi X/D near the center of the mask can scatter enough light into the field of view 
to scuttle a planet search (Kuchner & Traub 2002). Painting a graded-transmissivity mask 
requires a steady hand! This requirement has cast the classical coronagraph in an unfavorable 
light, despite its potential high performance and flexibility. 

In this paper, we offer a way around this pitfall of classical coronagraphy: an easy- 
to-manufacture class of image masks. We illustrate a family of binary image masks which 
offer a savings in construction tolerances of ~ 5 orders of magnitude compared to graded 
image masks, analogous to the advantage of using binary rather than graded pupil masks 
(Spergel 2001). These "notch filter" masks offer the same planet search and characterization 
advantages as ideal band- limited masks, providing a robust, practical means of controlling 
diffracted light in a planet-finding coronagraph. 



2. BAND-LIMITED MASKS 

We begin by reviewing the theory of band-limited image masks. We retain the notation 
of Kuchner & Traub (2002); image plane quantities have hats and pupil plane quantities do 
not. 

Figure 1 illustrates schematically how a coronagraph works; light passes through the 
pupil and converges on an image mask, then the pupil is re-imaged onto a Lyot stop. Starlight 
focused on the center of the image mask diffracts to the pupil edges, where the Lyot stop can 
block it, as shown on the left of the figure. Light from an off- axis planet diffracts all around 
the second pupil plane, as shown on the right of the figure, and largely passes through the 
Lyot stop. 

A band-limited mask has a transmission function chosen to diffract all the light from 
an on-axis source to angles within eD/{2\) of the edges of the pupil, as shown in Figure 1, 
so that a well-chosen Lyot stop can block identically all of that diffracted light. Such an 
image mask typically consists of a series of dark rings or stripes. The parameter, e, is the 
bandwidth of the mask. 

A mask can be described by an amplitude transmissivity, M{x,y), and intensity trans- 
missivity |M(a;, |/)p where x and y are cartesian coordinates in the image plane. Image 
masks are generally opaque (M = 0) in the center {x — y — 0) and close to transparent 
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Fig. 1. — Cartoon of a coronagraph with a band-limited image mask. The image mask 
diffracts on-axis starhght to a region restricted to the edges of the pupil plane, where a Lyot 
stop blocks it. Off-axis light from a planet diffracts all around the pupil plane, and through 
the center of the Lyot stop. 
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(M 1) away from the center, in the search area. To understand the need for band-hmited 
pupil masks, we must examine the Fourier transform of M{x, y), given by 

M{u,v)^ y"y"M(x,y)e-2"("^+''^)cixciy (1) 

The amphtude transmissivity of a completely transparent mask has only one Fourier com- 
ponent, at zero frequency, i.e. M{u,v) = 6{u,v). 

Figure 2 illustrates the operation of a mask with one cosine component besides the zero- 
frequency component, the sin^ mask (sin^ intensity transmissivity) described in Kuchner & 
Traub (2002). The Fourier transform of the amplitude transmissivity of this mask consists 
of three delta functions: 

M{u, v) = ~S{u + eD/{2X),v) + -S{u, v) - ^S{u - eD/{2X),v) (2) 
4 2 4 

This mask is the simplest example of a band-limited mask. 

The amplitude transmissivity, M{x, y), multiplies the field amplitude in the image plane. 
In the pupil plane, on the other side of a Fourier transform, this multiplication becomes a 
convolution. Figure 2a illustrates the convolution of the amplitude of the pupil field of an 
monochromatic on-axis source and the function, M{u,v), given in Equation 2. 

In the convolution, each 5-function from Equation 2 generates a weighted copy of the 
pupil field — a virtual pupil. We represent each copy of the field as a circle filled with + signs 
or — signs. The circular shape represents a circular aperture, though any aperture shape 
will do. Since the central 5-function has twice the weighting of the other 5-functions, the -|- 
signs have twice the density of the — signs in Figure 2a. 

Figure 2b depicts the field in the second pupil plane, the sum of the three virtual pupil 
fields shown in Figure 2a. In the center of Figure 2b, the densities of -|- and — signs are 
equal; for every + sign, there is a — sign. In this region, the fields cancel to zero. Elsewhere 
the fields do not cancel. The next optical element in the coronagraph beam train is a Lyot 
stop, which transmits light in the center of the pupil plane, but blocks the regions where the 
fields do not cancel. 

A given Lyot stop blocks the hght diffracted by a range of Fourier components. If the 
the Lyot stop blocks a fraction, e, of the pupil radius at the pupil edges, it will block the 
diffracted light from all spatial frequency components in the mask with spatial frequency 
\u\ < eD/{2X), where D is the telescope diameter, and A is the wavelength. One can create 
a mask which contains any or all of the cosine Fourier components at these low frequencies 
which the Lyot stop will still match; this family of masks which has power in only a limited 
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±+±+±±±±±±±±±±+±+± 
±+±+±±±±±±±±±±+±+± 
±+±+±±±±±±±±±±+±+± 
±+±+±±±±±±±±±±+±+± 
±+±+±±±±±±±±±±+±+± 
±+±+±±±±±±+±+± 
++++++++++ 



Fig. 2. — Convolving the pupil field of an on-axis source with the Fourier transform of sin^. 
a) For each delta function in the Fourier transform, there is one weighted copy of the field, 
a virtual pupil, b) The field in the second pupil plane is the sum of the fields of the virtual 
pupils. The fields cancel to zero in the center where there are equal densities of plus and 
minus signs. 
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range of low spatial frequencies is the family of band-limited masks. We can use e to refer 
to the bandwidth of a Lyot stop or the bandwidth of an image mask matched to that Lyot 
stop. 

Likewise, a given ideal Lyot stop and mask combination can work at a range of wave- 
lengths. The bandwidth of a given image mask is proportional to A, but the bandwidth of a 
given Lyot stop is independent of A. Therefore, a given Lyot stop/image mask combination 
will work at all wavelengths shorter than the wavelength for which it was designed. However, 
it can only have optimum throughput at one wavelength. 

Kuchner h Traub (2002) display a variety of one-dimensional band-limited mask am- 
plitude transmissivity functions. A useful compromise between search area and throughput 
is M(x) = (1 — sinc(7r,xeZ}/A)), where sinca: = s\\i[x)/x, and 1 — is the minimum 
value of sine a; (i.e., = 0.8215349763788L..). The throughput of a Lyot stop matched to a 
one-dimensional mask function is roughly 1 — e. Band-limited masks with additional Fourier 
components in the v direction are also possible, though to use these masks, one must stop 
the top and bottom of the pupil plane as well as the left and right. The throughput of such 
a Lyot stop is roughly (1 — e)^. 

At the request of NASA, a university-industry team associated with Ball Aerospace 
and Technologies Corporation studied a design for a space-based visible-hght planet finding 
telescope using a single 4 m by 10 m elliptical primary mirror. This team estimated that 
with a classical coronagraph using a Gaussian image mask, the design could detect an Earth 
twin orbiting a G2 V star at a distance of 10 pc in 0.86 hours, including time for 2 rotations 
of the image plane (Beichman et al. 2002). Once the location of the planet was known, a 
water band in the planet's atmosphere could be detected spcctroscopically in 0.14 days, and 
an O2 band could be detected in 0.8 days. If a single band-limited image mask of the form 
M(x) = 1— sinc^x were used instead of a Gaussian mask, the Lyot stop could be substantially 
widened, increasing the throughput, and the detection and characterization times would be 
reduced by a factor of roughly 0.7, to 0.6 hours for detection, 0.1 days for H2O, and 0.6 days 
for 02- 

Ideally, a band-limited mask combined with a Lyot stop completely blocks all on-axis 
starlight from reaching the second image plane, and attenuates off axis starlight to an easily 
manageable level. However, Kuchner & Traub (2002) discuss two significant limitations on 
the band-limited mask performance: pointing errors and errors in mask construction. None 
of the time estimates in the Ball report accounts for either of these errors, which affect all 
masks, band-limited or not. 

A mildly apodized Lyot stop can compensate for the leakage due to pointing errors 
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(Kuchner & Traub 2002). Apodizing the Lyot stop carries a throughput penalty, but even 
with this loss, the ideal classical coronagraph outfitted with a choice of band-limited masks 
remains by far the fastest of the idealized designs described in the Ball report for planet 
detection. Mask errors are more serious; all graded image masks suffer from impractically 
tight construction tolerances. We will show how to dramatically loosen the construction 
requirements by building binary masks. 



3. NOTCH FILTER MASKS 

To build a binary mask that retains the advantages of band-limited masks we will need 
to use more of the available function space for mask design. Section 2 reviewed the utility of 
masks whose Fourier components are limited to spatial frequencies < eD/{2\). However, 
there is another range of spatial frequencies available for mask design: as Kuchner & Traub 
(2002) described in their discussion of mask errors, high spatial frequency terms that diffract 
light well outside the opening in the Lyot stop do not affect the performance of a mask. We 
can add high spatial-frequency terms, with \u\ > (1 — e/2)(D/A), to the mask amplitude 
transmissivity function without altering the light admitted through the coronagraph as long 
as 

reD/{2\) 

/ M{u) du = 0. (3) 

J-€D/{2X) 

Figure 3 shows that the spatial frequency response of a general image mask which can 
completely suppress on-axis light resembles the spectral response of a notch filter. 

We can use the degrees of freedom available at high spatial frequencies to design masks 
which are relatively easy to construct to the necessary tolerances. For example, the transmis- 
sivity of a band-limited mask is analytic, so it can not have a constant value over any finite 
region. However, the transmissivity of a notch filter mask need not obey this restriction. 
The remainder of this paper will be a discussion of notch filter masks that take advantage 
of this opportunity. 



4. ONE-DIMENSIONAL MASK FUNCTIONS 

Sampling a function forces its Fourier transform to be periodic. We can harness this 
aliasing effect to generate useful notch filter mask functions. We will illustrate this principle 
first by considering mask functions of one variable only. Such a mask function can be realized 
as a striped mask as shown in Kuchner & Traub (2002). These functions can also be used 
as parts to construct two-dimensional masks. 
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Spatial Frequency 



Fig. 3. — Power spectra of a band-limited mask (a) and a notch filter mask (b). The mask 
functions may have power at spatial frequencies indicated by the hatched regions. 
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Throughout our discussion, Mbl{x) will be a function which can serve as the amplitude 

transmissivity of a band-limited mask: < Mbl{x) < 1, Mbl(O) = 0, and the Fourier 
transform, Mbl{u), of this function only has power at spatial frequencies u < eD/{2X). 
Such a function automatically satisfies Equation 3. We shall use Mbl{x) to create notch 
filter functions, Mnotch{x), that mimic Mbl{x) at low spatial frequencies. 



4.1. Sampling 

Actual image masks are constructed using finite-sized tools offering limited contrast and 
spatial resolution. We can design notch filter masks with this contraint in mind. Multiply 
Mbl{x) by a comb filter with spacing Ax to get a sampled version of Mbl{x), and convolve 
the result with a function P{x), to get 



Msampiedix) = P{x) * (^Mbl{x) Ax ^ 5(x - (n + C) Ax) ^ (4a) 
Msampiediu) = P(k) |^Mbl (^i) * ^ 5(k - n/ Ax)e-2-<^- ^ , (4b) 

where n ranges over all integers and * indicates convolution. For generality, we have allowed 
the sampling points to be offset from the mask center by a fraction ( of Ax. The kernel, 
P{x), can represent the "beam" of a nanofabrication tool. It should be normalized so that 
J^^P{x) dx = 1, and P{x) must be everywhere < l/(Aa;), so Msampied{x) remains < 1. 

This function we have created, Msampied{U') , only has power at \u — n/Ax\ < e/2. Its 
power spectrum resembles Figure 3 as long as the spacing between the samples satisfies the 
requirement 

Ax < \/D. (5) 

If C 7^ 0, then Msampied{u) typically has an imaginary component. However, if Equation 5 
holds, Mgampiediu) is always purely real at low frequencies {u < eD/{2\)). 

In general, Msampiediu) docs not match Mbl{u) at low frequencies, because the envelope 
function, P{u), in Equation 4b does not generally equal unity over the whole bandwidth of 
Mbl{x)- Rather, the envelope function tends to cause Msampied{u) to violate Equation 3. 
However, we can often correct for this effect and create a function, Mnotchi'^)-, which satisfies 
Equation 3 by subtracting a constant, Mq, from Msampiedi^). I.e., 

Mnotch — M sampled — Mq, (6) 
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where 

/eD/{2\) j-oo 
Msampled{u)du= / Mbl{ u)P{u) du — Mbl{x) * ■ (7) 

■eD/(2A) J-oo ^=0 

To use this technique, we must not sample Mbl{x) where Mbl{x) = 0, or else we 
will end up specifying a mask with Mnotch{x) < 0. Specifically, \(\ must be greater than 
some minimum value, Co, given by the condition that Mbl{Co^^) — -^o- We may sym- 
metrize the mask if wc desire by forming a combination such as {Mnotch{x) + M'^^^^f^{x))/2, 
or Mnotch{x)Mnotch{^) whcrc for Mnotchi^)^ ~C been substituted for (. The latter com- 
bination has twice the bandwidth of the former. 

For example, if we choose Ax — \/D, and P{x) — {D / \)Il{xD / \) , where n(x) is a 
tophat function, 

nW = li f (8) 

elsewhere, 

then P{u) = sinc^TiuX / D) , and our sampling algorithm generates a mask resembling a 
histogram plot. A graded version of this mask would consist of stripes of different uniform 
shadings. Choosing C = Co will generate a striped mask whose darkest stripe is perfectly 
opaque. Choosing Co > C ^ 0.5 will generate a striped mask which never becomes perfectly 
opaque, a potentially useful trick since fabricating graded masks with high optical densities 
can be a challenge (Wilson et al. 2002). 

If Mbl{x) = sin^ {nxeD/ {2 Xj), then Equation 7 tells us that for this mask, Mq = 
(1/2)(1 — sinc(7re/2)), and Co is given by the condition Mbl{CoX/D) = Mq = sin2(7rCoe/2). 
The trick probably only works for the sin^ mask when e — 1/n, since sin^ has so many zeros. 
Table 1 hsts Mq and Co for several other masks, given a tophat kernel. 

Figure 4 shows five examples of notch filter masks all of which are different versions of 
the same basic 1 — sine mask. Figure 4a is a simple band-limited mask with no additional 
high-frequency components. Figure 4b shows a version of this mask sampled as described 
above with a tophat kernel of width X/D. This kernel is the narrowest one that works with 
Ax — X/D; narrower kernels require finer samplings. 



5. BINARY MASKS 



In two dimensions, we can use the additional degrees of freedom afforded by the high- 
frequency terms in a notch filter function to generate a completely binary mask, i.e., a mask 




Fig. 4. — Close-up view of the center of a one-dimensional linear graded band-limited mask 
(a), and four notch-filter versions of this design (b,c, d, and e). 
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Table 1: Sampled Mask Parameters for P{x) = [D / X)Yi{xD / X) 



Mbl{x) 


Mo(e) 






e 


Mo 


Co 


^• 2 TTxeD 
2A 




sinc(7re/2)] 




0.2 
0.4 


0.0081842 
0.0322554 


0.2883579 
0.2873981 


N [1 sinc'^^;^^] 


N 
2 


1 - 


- 1 JJ^^ sinc(7ru') du' 




0.2 
0.4 


0.0022456 
0.0089032 


0.2039059 
0.2032511 


1 - (sinc^)^ 


1 
2 


1 - 


ii;^'{l-'-^)suic{nu')du'\ 


0.2 
0.4 


0.0013681 
0.0054400 


0.1019893 
0.1017713 



N < 0.82153497637881... 



which everywhere satisfies Mf,inary{x,y) = or Mbinaryix,y) — 1- Such a mask can be 
constructed entirely out of material which is highly opaque, like metal foil. 



Let 



5.1. Linear Binary Masks 



strzpeK ,yj >, ^ elscwhcre ^ ' 



and 



Mbinary{x,y) = ^^2^ 6 {y - uX / D) j * M stripe{x , y) . (10) 

The Fourier transform of this convolution is a product: 

Mbinaryiu, v) ^ (^J^ " uD / X)^ J J Mstripe{x, y)e-^-^^--+-y) dx dy 

= i^5{v-nD/X)^ J sinc{nv{X/D)M„,t,hix))Mnotch{x)e-''^'''''dx. (11) 
At low and mid-spatial frequencies, only the v — term contributes, and 

Mbinary{U: v) = Mnotchiu) ^ Mbl{u) for U,V < cD / (2A). 

Though Mnotch{x) = 1 in some places, it is possible to multiply Mnotch{x) by a positive 
real constant, less than 1, to allow for a mask substrate that is not perfectly transparent or 
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reflective, or to guarantee that the metal strips maintain at least a minimum width, at a 
small cost in throughput. 

If we use a sampled mask function for Mnotchi^), the binary mask may be constructed 
entirely from opaque rectangles of varying lengths, for example, generating a "manhattan" 
pattern for simple nanofabrication. Figures 4c, d and e show examples of binary masks which 
mimic the 1 — sine mask. Figures 4d and e are binary sampled masks. 

5.2. Circular Binary Masks 

We recommend using a linear mask for the following reasons: 1) Linear masks have 
bandwidth in one direction only, so they generally have the best throughput. 2) If one 
region of the mask deteriorates, the mask may simply be translated so that the starlight falls 
on a new region. 3) Errors in the uniformity of the wavefront in the direction perpendicular 
to the image mask cancel out in the Lyot plane; for example, the telescope need only be 
pointed accurately in one direction. 4) It may be possible to use a carefully oriented linear 
mask to block the light from a binary star. 

However, circularly symmetric image masks can provide slightly more search area than 
linear image masks, so we discuss them here. Figure 5a shows the center of a graded band- 
limited mask of the form MBiif^) = N{1 — sine irreD/X). Figure 5b shows a sampled version 
of this mask, where, necessarily, the sampling has been performed in two-dimensions. Cre- 
ating this sampled mask requires following the same procedure illustrated in Section 4.1 to 
guarantee that the mask satisfies Equation 3. The sample points are shifted by a fraction of 
X/D in some direction, and a constant is subtracted from the mask amplitude transmissivity. 



a) 




Fig. 5. — Close-up view of the center of a one-dimensional radial band-limited mask (a), and 
two equivalent notch- filter versions (b and c). 
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We can also replace an azimuthally symmetric transmission function, M{r), with a 
discrete K-fold symmetric "star" mask. First, choose a 1-dimensional band- limited function, 
M{r)BL or a notch filter version, Mnotchif)- Then let 

M (r- ft\ where mod(i^^/27r, 1) < Mnoteft(r) , . 

M,_,(r,e) = | ^ elsewhere ^^^^ 

The Fourier transform of this function (see, e.g., Jackson (1975), pl31 [problem 3.14]) is 

/•oo 

Mbinary{q,(p) = M{r)Jo{qr)r dr 
Jo 

m=o o „oo 

+ i"^ exp{imK(l)) / Jjnxiqr) sm{mMnotch{r))r dr (13) 
, ^ Jo 



m=— oo 



where q and (f) arc the radial and angular polar coordinates in the pupil plane, and J„ is the 
Bessel function of order m. Figure 5c shows an example of such a binary star mask. 

For a truely band-limited mask, the radial integrals in Equation 13 should be evaluated 
over a range from to infinity. However, as Kuchner & Traub (2002) discussed, a mask 
truncated at a radius of say, r = lOOA/D, can serve more than adequately as an approxi- 
mation to a band-limited mask. Moreover, the mask illumination falls off rapidly with r, so 
deviations from an ideal mask are inconsequential outside some radius rmax, which is likely 
to be much less than lOOX/D. 

If we consider the mask to be truncated at r = r^ax, the high frequency terms are 
significant only for high q. Since JmK{q^) ~ (qr)"^^ /{^"^^ {mK)\) for qr < mK, the higher 
order terms must have absolute values less than (Drmax/{4:X))^ /K\ inside the pupil {q < 
D/{2\)); with enough points in the star, these terms are all small. For example, if r^ax — 
lOX/D, M will be less than 10~^ interior to the Lyot stop for K > 14; this level suffices to 
allow a coronagraph to suppress the intensity of an on- axis source in the final image plane 
by a factor of 10~^°. 



5.3. Combining Notch Filter Masks 

In general, the product of two notch filter mask functions is not a notch filter mask 
function. However, all of the examples of notch filter mask functions discussed in this 
paper — except for the circular masks — have periodic Fourier transforms. The product of 
two such functions is another periodic notch filter function. For example, one notch filter 
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mask is 

f where \y-n\/D\ < M^otchix)\/i2D) 

Mbinaryix, y) = I OI \x - UlX/ D\ < Mr,otch{x)X/ {2D) (14) 

I 1 elsewhere. 

In such a product, the band widths of the component masks add in each direction separately. 

One may also produce a notch filter mask by combining the complements of periodic 
notch filter masks. For example, 

f 1 where \y - nX/D\ > (1 - Mnotch{x))X/{2D) 
Munary{x,y) =1 OY \x - mX/ D\ > (1 - Mnotch{x))X/ {2D) (15) 

I elsewhere. 

Figure 6 shows a close up of a mask with amplitude transmissivity (1 — sinc^a;)(l — sinewy), 
and a binary notch-filter version of this mask created by combining the complements of two 
mask functions of the form (1 — sinc^). This band-limited mask has a search area which 
closely resembles that of a common mask which is not band-limited — the Gaussian spot. 
This mask and the one in Equation 14 have bandwidth in both the x and y directions. 

As a third example, we can combine masks with sampled versions of the uniform mask, 
M{x) = C, a constant. Let 

Mt{x) = P{x) * Ax ^ 5{x - nAx) (16a) 

n 

Mt{u) = P(ii)^(5(ii-n/Ax), (16b) 

n 

where P{x) — C. We can multiply a notch filter mask function by Mt{x) and obtain 
another notch filter mask. If we choose Ax = X/D, and P{x) — {D / X)Il{xD / {C X) , then 
the new mask will look just like the old one, only painted with black stripes of width CX/D, 
spaced hj X/ D (Figure 4c), which may run in any direction. Since Mt{x) diffracts some light 
outside the Lyot stop, the intensity throughput of this new Mnotch{x) will be reduced — by 
a factor of C^. Combining binary masks and these striped masks may make it possible to 
design a range of masks which require no supportive substrate. 



6. MASK ERRORS 



Consider a binary mask like the one shown in Figure 4d, constructed from rectangles of 
opaque material, of width X/ D. How sensitive is the coronagraph to errors in the construction 
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of this mask? What if one of these rectangles, in the center of the mask, were accidentally 
made too short by an amount hX/D, where h « 11 

A missing rectangle of material — or an extra rectangle of material — would act like a 
tophat mask, diffracting light around the second pupil plane. A tophat mask is not band- 
limited and it has a power spectrum that falls off quite slowly with spatial frequency. There- 
fore it affords only modest cancellation of hght in the center of the second pupil plane. 

A tophat mask of width X/ D and length hX/D produces a diffraction pattern with most 
of its power in a zone with dimensions D/Xhy D/{hX). The intensity in this illuminated 
region is proportional to h, but only a fraction ~ /i of the illuminated region falls in the 
center of the Lyot stop. In this portion of the illuminated region which falls in the center 
of the Lyot stop, the field is roughly uniform, but attenuated by roughly a factor of 2; the 
intensity is attenuated by a factor of four. Therefore, the final image will acquire an extra 
image of a point source in the center, with fractional intensity ~ 0.25/i^. 

We can easily tolerate leakages of ~ 10^^'^ of the starlight falling in the center of the 
final image plane. If we are to avoid leaks of greater than this magnitude, we must avoid 
making the rectangles too short, unless 0.25/i^ < 10^^'^, or /i < 1/3000. For a telescope 
with focal ratio /, this tolerance becomes A//3000, or typically ~ 0.02 /im, for A = 1.0 /xm, 
/ = 60. 

If the error is not in the center of the mask, but in the search area, we can tolerate less 
leakage intensity, but the light falling on the hole will be diminished in intensity by a similar 
amount, so the requirement on the size of the hole remains about the same. A hole far from 
the center, outside the search area, say at lOOX/D, need only have h < a. few percent, since 
the wings of the stellar image that fall on it are typically four or more orders of magnitude 
weaker in intensity than the core of the stellar image. 

The tolerances for binary mask construction given here fall within the reach of standard 
nanofabrication techniques. Mask defects acquired during a mission may yield to compensa- 
tion by the active optics a planet-finding coronagraph will require to correct wavefront errors 
throughout the system. If the hole in the mask described above had dimensions X/Dhy X/D, 
but had an intensity depth of only /, then the hole would cause a fractional stellar leakage 
of ~ 0.25/, as opposed to ~ 0.25h'^; the shapes of binary masks are much less sensitive to 
errors than the intensity transmissivities of graded masks. 
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7. 



A WORKED EXAMPLE 



To further illustrate the use of a notch-filter mask, let us design some notch filter masks 
for a circular 4m TPF coronagraph (see, e.g., Brown et al. 2002). We will assume a bandpass 
from Xmin = 0.66 /xm to Xmax = 1-0 yum, and a mask construction tolerance of 20 nm. Rather 
than describing the optics in terms of the dimensionless diffraction scale, X/D, we will use 
the phisical size of the diffraction scale in the focal plane, A/, where / is the focal ratio of the 
telescope. The search problem and the characterization problem call for different speciahzed 
image masks, based on different band- limited functions. We will design a search mask — a 
linear mask rather than a circular mask for the reasons enumerated in Section 5.2. 

If we choose a Lyot stop that works at Xmax, the coronagraph will provide equal or 
better contrast over the whole band. The half power point of the mask, 9i/2, is angle from 

the optical axis where ^Mbl(^i/2)^ = 0.5; we will choose 9i/2 — ?>Xmax/D or 150 mas. 
This mask will enable us to search for planets as close as ~ 1 AU projected distance from a 
star at 6.5 parsecs. Searching for planets calls for a mask based on a band-limited function 

with small wings, providing a large search area where the planet image will be relatively 
unattenuated by the mask. A suitable one with the prescribed half-power point has the form 
Mbl = 1 - sinc^(7ra;e/(A„a^/)) with e = 0.4 at A„,„^, i.e., M = 1 - sinc^(a;/(1.76A„a^/)). 
Since the primary is circular, the shape of the Lyot stop will be the overlap region of two 
unit circles whose centers are separated by e, as depicted in Figure 4d of Kuchner & Traub 



To estimate the stellar leak due to pointing error and the finite size of the stellar disk, 
we may re-write Equation 16 in Kuchner & Traub (2002) so that it applies to any linear 
mask which is roughly quartic inside its half-power point. The fraction of the starlight that 
leaks through the mask, L, is 



where 6^, is the angular diameter of the star, is the pointing error in the x direction, and 
is the half power point of the mask. This equation also applies to nulling interferometers 
with quartic nulls. However, the half power point of a nulling interferometer's fringe depends 
on wavelength. 

When the pointing error is somewhat larger than the typical stellar diameter (^^ = 
1.43 mas for a G star at 6.5 pc), the pointing error dominates the leak. If the pointing 
errors, A^, are distributed in a Gaussian distribution about A^ — 0, the mean pointing- 



(2002). 



L 



el + AS{Mfel + 128(A^)4 



(17) 
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error-dominated leak is ^ 

where CAe is the standard deviation of the distribution. If we require a mean leakage of 
~ 3 X 10^^ of the starlight falling in the center of the final image, the pointing error tolerance 
becomes (XAe < ^^1/2/ 100, or 1.5 mas. 

This leakage due to pointing error can easily be suppressed to the 10~^° level in the 
search area given a weakly apodized Lyot stop. For example, an apodization function of the 
form cos7rMA/((l — e/2)D) provides 2.5 orders of magnitude of suppression at 3X/D. The 
total throughput of the coronagraph would be 1 — e = 0.6 without the Lyot stop apodization; 
with the apodization realized as a binary mask, it is 0.358. This apodization is workable, 
but not optimal; further work on choosing matched pupil and image stops could improve the 
overall system performance. 

Clearly, the leakage due to pointing errors quickly shrinks for planet-finding corona- 
graphs with larger inner working distances, like the proposed Eclipse mission (Trauger et al. 
2002a). For this 2 m class telescope with inner working distance ~ 300 mas, a pointing error 
of (J AS = 3 mas suffices to match the above performance. Likewise, if the pointing errors and 
other low-spatial-frequency errors could be controlled below the levels assumed here, these 
coronagraph designs could operate at smaller inner working distances. 

We will realize the mask as a binary notch filter mask like the one in Figure 4d. In 

Section 6, we showed that the lengths of the bars in this mask must be accurate to A//3000. 
Since a hole of a given physical size subtends a larger fraction of X/D at smaller A, this 
tolerance applies at Xmin- In other words, if we can manufacture a mask accurate to 20 nm, 
we require a focal ratio / = 90. The strips must have maximum width < Xminf = 59.4/xm. 
Explicitly, the mask function would be: 

M,,^ary{x,y) - ^ 1 elscwhere. 



where 



Table 1 shows that our mask has Mq = 0.0054400 and C = 0.1017713. 
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In the absence of noise, analyzing the data from this coronagraph is trivial. A planet at 
angle 9p from the optical axis is simply attenuated by the low spatial frequency components 
of the notch ffiter mask, i.e., f=i (l — sinc^(^p/90mas))^. The shape of the PSF is set entirely 
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by the Lyot stop; it is the squared absolute value of Fourier conjugate of the Lyot stop 
amplitude transmissivity. 

8. SIMULATION OF CORONAGRAPH PERFORMANCE 

We numerically simulated the performance of this notch-filter coronagraph design by 
following the Fraunhofer propagation of light through a coronagraph using fast Fourier trans- 
forms. Nakajima (1994), Stahl & Sandler (1995) and Sivaramakrishnan et al. (2001) have 
used this technique to model the performance of ground-based coronagraphs. We simulated 
the broad-band performance by running the monochromatic simulation 10 times over a range 
of wavelengths from 0.66 to 1.0 fim and averaged the output images weighted by the stellar 
flux, assuming a Raylcigh- Jeans law star and planet. The noise representations remained 
the same from wavelength to wavelength — scaled appropriately to model pathlength errors 
rather than phase errors and to reflect the change in the diffraction scale. 

We used a 1024 by 1024 grid with resolution Amm/(2-D). In this representation, errors 
in the shape of the image mask become variations in the mask amplitude transmissivity. For 
example, if a bar in the binary mask were too long by hX/D, the values of four adjacent 
elements in the mask amplitude transmissivity matrix would increase by h. 

We assumed a circular pupil and a Lyot stop with the shape of two overlapping circles 
as described above multiplied by costtm/ ((1 — e/2)D)). Seen from afar, the contrast between 
the Earth and the Sun when the Earth is at maximum angular separation is 2 x 10~^° (Des 
Marais et al. 2001). We used this contrast level for the planets in our simulation. 

Figure 7 shows the intensity in the second pupil plane before (a) and after (b) the Lyot 
stop for a monochromatic simulation at 1.0 iim. of a system containing a star and a planet 
at 2QXjnax/D. The intensity after the Lyot stop has been multiplied by a factor of 10^. This 
figure represents a more accurate version of the cartoon in Figure 1. Figure 7a shows how 
the diffracted light falls within regions of width e around the left and right edges of the Lyot 
stop. The notch filter mask adds further illumination to the Lyot stop, farther off axis. Our 
simulation does not model these artifacts; the Lyot stop blocks them completely. The planet 
adds a uniform background illumination to the region inside the Lyot stop, though the noise 
peaks due to wavefront errors and mask errors dominate Figure 7b. 

Figure 8 shows a cut through the final image plane along the x-axis. The top panel of 
the figure shows a plot of log (l — sinc^(6'p/90mas)) . This quantity is the intensity attenu- 
ation supplied by the coronagraph, neglecting any modification of the low spatial frequency 
components of the mask that might occur in the construction of a notch filter representation. 
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Fig. 6. — Close-up view of the center of a (1 — sinc^a;)(l — sinewy) mask, a band-limited mask 
which offers nearly as much search area as a Gaussian spot, (a) Graded version (b) Binary 
notch filter version. 



«) 
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Fig. 7. — Simulated monochromatic intensity in the Lyot plane for a star and planet system 
imaged by the coronagraph described in Section 7 a) before the Lyot stop b) after the Lyot 
stop, amplified by a factor of 10^. 
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The lower panel of the figure shows the relative surface brightness of several noise 
contributions to the image, normalized to the peak of what the stellar image would be if 
the image mask were removed. The dotted curve shows the 10~^'^ leakage due to pointing 
error which we prescribed. The form of this curve is simply the Fourier transform of the 
Lyot stop averaged over the band weighted by the stellar flux. The dash-dot curve shows the 
consequence of adding white noise to the lengths of the bars with r.m.s. 20 nm. This noise 
concentrates near the optical axis in the final image plane because it is multiplied by the 
image of the star that falls on the image mask. This phenomenon suggests that the tolerance 
of the coronagraph to mask errors can be altered, and possibly improved by manipulating 
the shape of the entrance pupil. 

The solid black curve shows the consequence of adding white-noise phase and amplitude 
errors to the incoming wavefront: fractional amplitude errors of r.m.s. 10~^ over spatial 
frequencies corresponding to the search area in the image plane (3-60 \/ D) and phase errors 
due to pathlength errors of r.m.s 0.5 A over these frequencies. A pathlength error of 0.5 A 
corresponds to an error in the flgurc of a mirror of 0.25 A. Preliminary tests in the High 
Contrast Imaging Testbed at JPL have demonstrated deformable mirror wavefront control 
at this level (Trauger et al. 2002a,b). Wavefront errors clearly dominate mask errors and 
pointing errors, except within a few diffraction widths of the optical axis. 

The grey curves show the images of two planets located at 3Xmax/D and 20\„iax/ D, 
i.e., 154 mas and 1026 mas — or 0.8 AU and 5.1 AU for a star 5 pc distant. The planet at 
^^max/D is attenuated by a factor of 0.5 because it sits at the mask's half power point. If Q 
is the contrast between the planet's peak intensity and the local scattered light background 
(Brown & Burrows 1990), the 3X/D planet has Q ^ 0.5, and the 20X/D planet has Q ~ 1. 
Planets with Q ~ 1 can easily be detected in a coronagraph using spectral deconvolution 
techniques, for example, given sufficiently low photon noise (Sparks & Ford 2002). 

9. CONCLUSION 

We have illustrated the use of notch filter functions to generate several kinds of image 
masks which should be relatively easy to manufacture. We showed graded masks whose 
transmissivities are everywhere greater than zero. We showed binary image masks, which 
can be cut or shaped from pieces of opaque material. These binary masks can be manufac- 
tured to the tolerances necessary for terrestrial planet finding using standard nanofabrication 
techniques, and can potentially be made self-supporting. Our simulations of the performance 
of a coronagraph outfitted with a binary notch filter mask suggest that this technique could 
reveal extrasolar planets similar in brightness to the Earth around nearby stars, given fore- 
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Fig. 8. — Broadband simulation of images produced by the notch filter coronagraph design 
described in the text. The dotted curve shows leakage due to pointing error. The dot-dash 
curve adds errors to the lengths of the mask bars. The solid curve adds amplitude and phase 
errors to the incoming wavefront. The grey curves show the images of planets with relative 
flux 2 X 10~^° at 3 Xmax/ D and 20 Xmax/ D. The upper panel shows the attenuation caused 
by the image mask. 
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seeable improvements in wavefront control on a highly stable space platform. 

Binary notch filter masks combine many of the advantages of binary pupil masks (ease 
of manufacture, achromaticity, robustness) with the advantages of band-limited image masks 
(large search area, and small inner working distance). Using binary pupil or image masks 
seems to inevitably require stacking many copies of the same basic aperture shape; Kasdin, 
Spergel & Littman (2001) used this principle to generate binary pupil masks; we have used 
it to generate binary image masks. In Kasdin, Spergel & Littman (2001), the high-spatial 
frequency artifacts of this stacking procedure appear in the image plane directed away from 
a search sector. In notch filter masks, the high-spatial frequency artifacts are directed into 
the Lyot stop. 

Ultimately, a space telescope for direct optical imaging of extrasolar planets may in- 
corporate more than one diffracted- light management strategy. Having a choice of different 
techniques available will allow a mission to adapt to changing observing needs as our un- 
derstanding of high-contrast space telescopes improves and the phenomenology of extrasolar 
planets unfolds. 
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